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Abstract 

We present a discussion on how to define the running electromagnetic cou- 
pling constant at Mz or some other intermediate scale as e.g. m~[. We argue 
that a natural definition consistent with general requirements of the renormal- 
ization group should be based on Euclidean values of the momentum of the 
photon propagator as the appropriate scale. We demonstrate in an explicit 
example of evaluating the running coupling constant at the scale of the T res- 
onance mass that the usual definition of the hadronic contribution with a prin- 
cipal value prescription is inconsistent. In the determination of the value of a 
at Mz the numerical difference due to using a Euclidean definition rather than 
the principal value one is comparable in size to the errors caused by existing 
experimental and QCD inputs to the evaluation of a(Mz)- 



In applications to high precision tests of the standard model [Q] with observables 
near the Z boson peak the electromagnetic coupling constant should be used at a scale 
of the order of the Z boson mass Mz, (see e.g. [|[ |J). The running electromagnetic 
coupling constant at Mz has even been chosen as a standard reference parameter [||. 
It differs numerically from the value of the fine structure constant a -1 = 137.036 . . . 
defined at zero momentum or from the Coulomb law for heavy nonrelativistic parti- 
cles. The change is usually accounted for through the renormalization group equation 
Because the fine structure constant is defined at vanishing momentum and is 
taken as initial value in the solution of the renormalization group equation, the run- 
ning electromagnetic coupling constant at Mz is an infrared sensitive quantity in 
as much as the contribution of strong interaction is not easy to compute due to 
the nonperturbative region at small energies. Therefore this contribution is usually 
taken into account in the leading order of electromagnetic interaction within a semi- 
phenomenological approximation through a dispersion relation. There has been a 
renewal of interest in a precise determination of the hadronic contribution during the 
last years in particular in connection with the constraints on the Higgs boson mass 
. Some recent references giving a state-of-the-art analysis of this contribution are 



10, 11]. A quasi-analytical approach was used in |T2j where some references 



to earlier paper can be found (see also |L3|, 0). An extremely thorough data-based 



analysis is given in 15 



In the present note we critically discuss the definition of the running electromag- 
netic coupling constant at Mz as it is used in the literature. The standard approach 
consists in using the principal value prescription at the appropriate scale in the phys- 
ical domain on the positive energy semiaxis. We argue that a natural definition 
consistent with general requirements and the standard notion of running used in 
renormalization group applications should be based on the Euclidean momentum of 



2 



the photon propagator as the appropriate scale. 

The running coupling a(q 2 ) is defined through the (one-photon irreducible) photon 
vacuum polarization function IT 7 (g 2 ) as 

Il 7 (g 2 ) contains both leptonic and hadronic contributions. The hadronic part of the 
polarization function Il 7 (g 2 ) (with one subtraction at zero momentum) reads 

OO 

n; a V) = -fg 2 / ( Rh{S 2 ds ^ ( 2 ) 

1 6Ti J s[s — q 2 — i0) 

where Rh{s) is the normalized cross section of e + e~ annihilation into hadrons. Let 
us introduce, for convenience, the polarization function n(g 2 ) 

ntf) = -e r . Ri T m p) 

Jo s(s — q 2 — i0) 



such that 



n 7 (g 2 ) = -I% 2 ) (4) 



where R(s) is the corresponding spectral density. Note that a(q 2 ) is defined for every 
complex value of q 2 by eq. (|J). For real negative q 2 the polarization function Il 7 (g 2 ) 
(and Il(g 2 ) as well) is a positive real number because the spectral density R(s) is 
positive. 

The definition (|I|) is used in renormalization group applications and the scale 
q 2 is taken to be a real negative number which corresponds to a propagator in the 
Euclidean domain. The Euclidean definition is mainly used in applications of grand 



unified theories ||16| , supersymmetry at large energy |L7|, physics at the Planck's scale 
etc. 

For the precise study of the physics at the Z boson pole the effective electromag- 
netic interaction coupling constant a is represented in the form 



Numerically one obtains a real positive number for Aa. 

A theoretical expression for Aa is defined in the present literature directly on the 
positive semiaxis by making use of the principal value prescription for the singularity 
of the integrand in eq. @ 



Here P / denotes the principal value of the integral. This makes Aa real (the initial 
iO prescription for the integral gives it an imaginary part) which is appropriate for a 
coupling constant. We argue that this prescription is not adequate for the physical 
situation at hand and does not correspond to the notion of a running coupling used 
in the standard renormalization group applications. The latter corresponds to scales 
taken in the Euclidean domain 



The running electromagnetic coupling constant at the scale Mz is defined then as 
otE{M z ). Therefore, for the phenomenological parameter "running electromagnetic 
coupling constant at the scale M z " denoted by a we have two representations: 
i) the standard one with a principal value prescription 



Aa = Ren 7 (M|), 










ii) the alternative one in the Euclidean domain 





(9) 



which is defined through 




(10) 
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We suggest that the Euclidean version is used. The idea of changing scales is 
embodied in the renormalization group equation which allows one to control large 
logarithms. Therefore theoretically one is dealing with a logarithm of the ratio of two 
scales. Note that the notion of scale becomes rather imprecise as soon as complex 
numbers are involved. For example, numerically M§ has the same 'scale' as e m M\ = 
— M§. The choice of an appropriate scale is determined by the particular kinematics 
and the higher order corrections of each particular process in question. In the leading 
logarithmic approximation, however, keeping the finite corrections to large logarithms 
is beyond the accuracy of the approximation that renders all scales with the same 
absolute value equivalent. As a reference value for the coupling constant the usual 
choice is to take a Euclidean point. This problem was discussed for strong interactions 
rj| where different versions of a real part or absolute value definition of the coupling 



in 



constant at complex points have been also studied. In the general case it is difficult 
to decide on how to deal with observables including complex numbers within the 
renormalization group resummation of logarithms. For two-point functions, however, 
there is a natural solution to this problem based on their analytic properties given by 
the dispersion representation |L9], pQ[ . 

Below we discuss these two possibilities of defining the "running electromagnetic 
coupling constant at the scale Mz" ■ 

First we show that the two definitions (Euclidean and principal value) are close 
numerically for applications in the vicinity of the Z boson peak discussed in the 
literature. Let us take eq. @ and split the whole region of integration into two parts 
separated by so 



W } q Jo s(s -q 2 - zO) q Jso s(s - q 2 - iO) 1 } 

If \q 2 \ is chosen such that \q 2 \ ^> Sq one can expand the denominator in the first 
integral. Then, if sq is large enough one can use perturbation theory for the spectral 



density in the second integral. For illustrative purposes we choose a very simplified 
approximation for R(s), namely R(s) = const = 1 for s > Sq. Then one obtains 



which is independent of the phase of the complex number q 2 . The same result can be 
obtained also directly from eq. (|]) in this limit. Therefore in the above approximation 
with the suggested regime of variables both the Euclidean and principal value defini- 
tions are equivalent numerically. Later on we discuss corrections to this leading order 
approximation which depend on whether the Euclidean and principal value definition 
is used. 

This is a qualitative picture. Because the above simplifying assumptions can be 
expected to correctly embody the main features of a more sophisticated numerical 
analysis it is clear that the numerical change stemming from the choice of q 2 in the 
Euclidean domain rather than using the conventional definition is under control at the 
scale of Mz and does not jeopardize current phenomenology. However, the definition 
in the Euclidean domain given in eq. (|H]) is preferable from a theoretical point of 
view. It is natural. It gives a real number. It is smooth. It is consistent with the 
renormalization group. 

And the principal value prescription has equally obvious deficiencies. It is ad hoc. 
It gives a real part of a propagator which is not directly related to a coupling constant 
in a renormalization group sense. It is not smooth. 

The last deficiency is, in fact, the most crucial one. Let us present more details. 
We take the principal value definition at q 2 = M\ and compute the polarization 




(12) 



where the principal value prescription has been used. Expanding eq. (|T2| ) in the limit 
\q 2 \ s one finally obtains 




(13) 
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function for a model spectral density R(s) = 9(s — so) obtaining 

II(Mf ) = In |S ° ~ M|1 . (14) 
So 

Note first that the polarization function (|I~4D gives a rather curious result 



OC 



PI 77-^ = (15) 



ds 

s(s ~MV) 

which means that the contribution of all states with masses larger than Mz/V2 ~ 
60 GeV is exactly equal to zero assuming the asymptotic spectral density in this 
region to be a constant. Also there is a sign change in the vicinity of Mz/V2- This 
feature persists for any realistic R{s) in the vicinity of some point s* ~ M§/2 because 
in this region QCD perturbation theory works well and the spectral density is rather 
smooth and close to its asymptotic value which is almost a constant (up to a slow 
logarithmic decrease). Therefore one gets the exact equality 

p/-*W* =0 (16) 
J s(s -Ml) y ' 

s* 

for some s* ~ Mf/2. 

Furthermore, if one takes So = M\ in eq. (|i4|) then the logarithm is ill-defined. 
This can be seen to be a consequence of the principal value definition eq. @. Even 
if this is a rather academic example we nevertheless take it as a warning (because 
there is no sharp increase of the spectral density or changes in general in the vicinity 
of the Z boson mass). More realistic situations are considered below. In contrast to 
the principal value definition the Euclidean definition is fine also in this case 

n(-Mf) = In f°±M . (17) 
So 



The reason for the ill-defined behavior of eqs. ( 14j) and (y) is clear. The principal value 



prescription leads to a distribution which is defined only on smooth functions. A 
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product of two distributions 



P 



s-M\ 



1 



) 



6(8 - Ml) 



(18) 



is not an integrable function. An ad hoc definition with a principal value prescription 
fails to define a value for the running coupling at some particular points and one has 
to introduce further rules for such cases. 

Also in a more realistic situation one needs the running electromagnetic coupling 
constant at the scales around the masses of resonances of the J/ip or T families to 
account for their leptonic widths. With the principal value definition it is impossible 
to compute the running electromagnetic coupling constant at the scales around the 
resonance masses. Indeed, the correction Aa to the running electromagnetic coupling 
constant at the scale of the T meson mass m T is given by an ill-defined integral of 
the product of two distributions 



This quantity is not defined as a distribution for the same reason as a product of two 
distributions is not defined when their singular points coincide. 

Leaving the mathematical statement about the ill-defined behavior of a product 
of two distributions aside, in practice (for finite widths of the resonances or for some 
sharp but still smooth increase of the spectral function) the results following from the 
principal value definition will be unstable. As an explicit example we take the spectral 
density corresponding a single Breit-Wigner resonance and calculate its contribution 
to iten(Mf). The Breit-Wigner spectral function is given by 



where M and T are the mass and the width of the resonance. As T — > one obtains 
Rbw(s) — > 5(s — M 2 ). After integration (with proper care for the point s = 0) one 




(19) 



Rbw(s) 



1 YM 



(20) 



7T (s - M 2 ) 2 + Y 2 M 2 
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finds 

M 2 p r RBw(s)ds _ M^-Ml _ A 

z J s{s- Mf) (M 2 - M|) 2 + T 2 M 2 A 2 + T 2 M 2 1 ; 

with A = M 2 — M\ . Only potentially singular terms for the limit r — > in the 
vicinity M 2 ~ Mf have been kept. The last function in eq. (|21|) has its extremal 
points at A = ±TM with the values 



A 



= ±— *— . (22) 
A =±rM 2fM v ' 



A 2 + r 2 M 2 

There is no regular limit r — > and further rules would be required to deal with 



this limit. Note that the result ([21]) can be obtained without explicit integration of 
the Breit-Wigner spectrum. Since the Breit-Wigner function can be regarded as a 
regularization of a 5-distribution one can introduce a regularization of P^ in the class 
of infinitely smooth function instead. An example is 

hm — = P- and hm — — nsr f ' - = P 



i-o x 2 + e 2 x e-o ( s - Mf) 2 + e 2 s - Mf ' 

Then after the integration with an infinitely narrow resonance one gets 

2 t 5(s - M 2 )( g - Mpds _ M 2 A A 

z 7 s[(a - Ml) 2 + e 2 ] M 2 A 2 + e 2 A 2 + e 2 1 J 

for M 2 ~ M\. The regularization can not be unambiguously removed, i.e. there is 
no unique limit at e = in the vicinity of A = 0. Of course, this is a reflection of the 
fact that the product of two distributions is ill-defined. 

We will not dwell on the ill-defined behavior when a ^-function type spectral 
density is used. This is a realistic situation when one computes the light quark 
contributions to the running electromagnetic coupling constant normalized in the 
vicinity of a sharp raise of the spectral density around 1.5 GeV 2 . 

With a Euclidean definition none of the above difficulties appear. Also because in 
this case the polarization function is defined in the Euclidean domain one need not 
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integrate all data (only a small region near the origin requires explicit integration) 
and can use all the power of perturbation theory (e.g. PI| ). Although special care 
has to be taken about the subtraction at zero momentum that enters the definition of 
the coupling constant and makes it an infrared sensitive quantity. For this purpose, 



however, more sophisticated means can be used that will increase the accuracy [22 



23, 24 



We add two further remarks. The first concerns the leptonic contribution. For a 
lepton with the mass mi the asymptotic form of this contribution reads 

'Mf\ 5 
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n ^ = ln Uf "3 (24) 



and has the same real part for any real phase <p of e lip M z . When the asymptotic form 
is used both prescriptions are equivalent numerically. 

The second remark is related to the higher order contributions of the electroweak 
interactions. In the next order of the electroweak interaction there is a contribution 
of the Z boson peak to the polarization function U^(q 2 ) due to transitions (e.g. 



25|j). Therefore in that order one has to interpret the product of the principal value 
distribution with the sharp Breit-Wigner spectrum of the Z boson pole itself. Again 
this problem is not present within the Euclidean definition. 

Now we discuss very briefly the numerical difference that can result from the 
change of the definition of the running coupling constant at Mz- Our model for 
the hadronic spectral density R(s) is simple and is mainly designed for illustrative 
purposes such that one can easily trace the difference between the principal value 
and Euclidean definitions of the running coupling constant. The spectral density is 
chosen such that all calculations can be done analytically which is convenient for the 
purpose of estimating the order of magnitude of the difference in the two definitions. 
For the light quarks u, d, s we assume the existence of a low lying resonance (like p, 
u, ip) and a continuum. In general, we take the following form of the spectrum for 
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every light quark flavor 

Rught(s) = 3Q 2 q [2m 2 5{s - m 2 ) + 6{s - 2m 2 



according to the model of ref. [26] with Q q being a light quark fractional charge. The 
coupling of the low lying resonances have been replaced by the duality interval 2m 2 . 
For heavy quarks we take the simplest model of the form 

Rheavy(s) = 3Qq6»(s - AuIq) 

which represents the partonic asymptotic value with a naive step function two-quark- 
threshold with Qq being a heavy quark charge. Collecting everything together we 
find the following results. The three light quarks give the result for II(g 2 ) in the 
general form 

/ 9x ( -2<7 2 , |2m 2 -g 2 |\ 

lWg 2 )=2( I ^ +ln l 



Expanding at \q \ ~ M§ one gets 



n 



light 



m 2 — q 2 — iO 2m 2 



2 + ^ M + 0(m 6 /Mf) 



(25) 



2m 2 

which gives the same answer for both definition with a high precision. The difference 
starts at the order 0(m 6 /M|) and is completely negligible for light resonance masses 
of order 1 GeV (p-, u-, (^-resonances, for instance). For the charm contribution we 
find 

PV A M 2 -Ami 

charm 3 ^ 

in the case of the principal value prescription. In the case of the Euclidean prescription 
one has 

4 Mf + 4mg 

charm 3 ^ ■ 

Expanding these formulae in the small ratio 4m 2 /M§ one finds in the leading order 



pv,e 4 M 2 1 



T r r V ,-C/ 1 

ll charm "3 4 m 2 



4m 2 1 

1 =F 



Mi ]n Mi 
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We keep this form for further numerical comparison in the case of b and t quarks. For 
the b quark contribution one gets 



1.4 GeV, m b = 4.8 GeV, m t = 175 GeV, M z = 91 GeV. Note that the exact 
definition of the quark mass parameters is not required here because it is far beyond 
the accuracy of our simple model. Even more, these parameters can be considered 
as effective parameters serving to describe integrals over the threshold regions of 
quark production. Nevertheless we stick to almost canonical values for the pole 
quark masses. While the absolute value of the contribution will be obtained rather 
approximately, the model however is fairly sufficient for our main purpose to estimate 
the difference between the two definitions. Numerically, eq. (|25| ) leads to the light 
quark contribution 




and for the t quark 




For numerical estimates we take \/2m = 1 GeV for the light quark resonances, m c 



light quarks (u , d , s) = 2 x (2 + 2 x 4.5) = 22.0 



with both prescriptions. For the contributions of the heavy quarks we find 



b — quark 



c — quark 



t — quark 



- x 7.0(1. =F 0.14 x 10~ 3 ) 
3 

- x 4.5(1. =p 2.5 x 10~ 3 ) 
3 

- x (^68 x 10~ 3 ). 



(26) 



Summing everything together one obtains 



4 14 
22.0 + - x 7.0(1. =f 0.14 x 10~ 3 ) c + - x 4.5(1. =f 2.5 x 10" 3 ) fe + - x (^68 x 10" 3 ) 

o o o 
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= 32.8 + 0.15^33 + 0.15 (27) 

where 5 = — 1 for the principal value definition and 5 = 1 for the Euclidean definition. 
One sees that the difference is saturated by the top quark contribution. This is natural 
because it has a mass closest to Mz- Its contribution is small in absolute value but 
is completely different for the two definitions. 

The lepton contribution is taken into account according to the asymptotic formula 
eq. (p4|) with m e = 0.5 MeV, m M = 0.1 GeV, m T = 1.8 GeV which gives 

(24.2) e + (13.6) M + (7.8) T - 5 = 40.6 « 41 . 

The final result for the total contribution of the charged fermions to the polarization 
function reads 

33 + 0.15 + 41 = 74 + 0.15. 
Dividing by 3tt one obtains the total contribution to the inverse coupling constant 

(74 + 0.15)/3tt = 7.9 + 0.015. 

For the inverse running electromagnetic coupling constant we obtain 

137.0 - (7.9 + 0.01 5) = 129.1 - 0.01 5 . 

Even though the central value of our approximate evaluation is rather close to the 
results of more precise evaluations || 10, |ll| this agreement should not be taken 



too seriously. Our estimate is rather rough and serves the purpose well to obtain 
the numerical change between the two definitions of the running electromagnetic 
coupling constant at Mz- At present the change due to the different definitions has 
no influence on current phenomenology. This change is within the error bars for the 
uncertainty of the more precise values 128.93 ± 0.06 |J or even with smaller errors 
128.93 ± 0.015 cxp ± 0.015th g 0, 0. Therefore the Euclidean definition, which 
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we consider more consistent theoretically, does not violate current phenomenology. 
However, the change is comparable in size with the present uncertainties and when the 
experimental data used in the determination of the running electromagnetic coupling 
constant at Mz improves in the future the difference between the two definitions will 
become significant. 

Our last numerical example concerns the order of magnitude of the singular term in 
the coupling normalized at the scale mx within the principal value prescription. With 
the same normalization as in our simple model the contribution of the T-resonance 

to the spectral density in the Breit-Wigner approximation reads 

2 

Ry{s) = -m r A r R BW (s,m r ,r) (28) 



where Ax ~ 1 GeV is its duality interval in energy units [f27| related to its leptonic 
decay width, A T = 277rr(T -> e + e-)/2a 2 with T(T -> e + e~) = 1.32 keV while 
T = 52.5 keV is its full width. The singular contribution of the T-resonance to the 



polarization function abruptly changes between the two extremes taken from eq. (g2| 

2 A / 1 
±-mxAx 



"3 V2rm T , 

as the normalization point passes the position of the resonance. Numerically one gets 

A T □ 
— = 6.7 x 10 3 

3r 

which is far too big from the point of view of phenomenology. 

The examples presented in this paper clearly demonstrate the inconsistency of 
the present definition of the running electromagnetic coupling constant within the 
principal value prescription. These problems are not noticeable, however, when one 
discusses a normalization point around Mz because the hadronic spectral density is 
smooth in this region. 

We are going to present the results of an accurate numerical analysis within the 
Euclidean definition elsewhere. 
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To conclude, we suggest that the Euclidean domain definition of the running 
electromagnetic coupling constant at Mz is used as a reference parameter for high 
precision tests of the standard model at the Z boson peak. It is free of the shortcom- 
ings of the present definition which is based on the propagator at the physical value 
of the Z boson mass within the principal value prescription. 
Acknowledgements 

The work is supported in part by Volkswagen Foundation under contract No. 1/73611. 
A. A. Pivovarov is supported in part by the Russian Fund for Basic Research under 
contracts Nos. 96-01-01860 and 97-02-17065. The present stay of A. A. Pivovarov in 
Mainz was made possible due to Alexander von Humboldt fellowship. 

References 

[1] Precision tests of the standard electroweak model. Ed. P. Langacker, World Sci- 
entific, 1995. 

[2] D. Bardin et al., Electroweak working group report. CERN-95-03A, Sep 1997. 
158pp. 

[3] W. Hollik, J. Phys. G23 (1997) 1503 

[4] Particle Data Group, Eur. Phys. J. C3 (1998) 1 

[5] E.C.G. Stueckelberg, A. Petermann, Helv. Phys. Acta 26 (1953) 499 

[6] M. Gell-Mann, F.E. Low, Phys. Rev. 95 (1954) 1300 

[7] J.Gunion et al., Frontiers in Physics, Addison- Wesley, Menlo Park, 1990 

[8] M. Davier, LAL-98-87, Dec 1998. 13pp. |hep-ph/9812370 



15 



[9] S. Groote, J.G. Korner, K. Schilcher, and N.F. Nasrallah, 
Phys. Lett. B440 (1998) 375 

[10] J.H. Kuhn, M. Steinhauser, Phys. Lett. B437 (1998) 425 

[11] M. Davier, A. Hocker, Phys. Lett. B435 (1998) 427 

[12] N.V. Krasnikov, R. Rodenberg, Nuovo Cim. 111A (1998) 217 

[13] V. Novikov, L. Okun, M. Vysotsky, A. Rozanov, 
ITEP-PH-6-97, Jul 1997. 13p, |hep-ph/9707359 

[14] H. Burkhardt and B. Pietrzyk, Phys. Lett. B356 (1995) 398 

[15] S. Eidelman and F. Jegerlehner, Z. Phys. C67 (1995) 585 

[16] P. Langacker, Phys. Rept. 72 (1981) 185 

[17] H.P Nilles, Phys. Rept. 110 (1984) 1 

[18] A. Bottino and C.W. Kim, Lett. Nuovo Cim. 30 (1981) 353; 
Nuovo Cim. 57A (1980) 459; 

N.V. Krasnikov and A.A. Pivovarov, Phys. Lett. 116 B (1982) 168; 

M.R. Pennington and G.G. Ross, Phys. Lett. 102 B (1981) 167; 

M.R. Pennington, R.G. Roberts and G.G. Ross, Nucl. Phys. B242 (1984) 69 

[19] A.A. Pivovarov, Sov. J. Nucl. Phys. 54 (1991) 676; 

Z. Phys. C53 (1992) 461; Nuovo Cim. 105A (1992) 813 

[20] S. Groote, J.G. Korner, A.A. Pivovarov, 

Phys. Lett. B407 (1997) 66; Mod. Phys. Lett. A13 (1998) 637 

[21] K.G. Chetyrkin, J.H. Kuhn and A. Kwiatkowski, Phys. Rep. 277 (1996) 189 



16 



[22] M.A. Shifman, A.I. Vainshtein and V.I. Zakharov, Nucl. Phys. B147 (1979) 385 

[23] N.F. Nasrallah, N.A.Papadopoulos and K. Schilcher, 
Z. Phys. C16 (1983) 323 

[24] K.G. Chetyrkin, A. Kwiatkowski, Z. Phys. C59 (1993) 525 

[25] See W.Hollik contribution in 0, p.37. 

[26] N.V. Krasnikov, A. A. Pivovarov, N. Tavkhelidze, 
Z. Phys. C19 (1983) 301 

[27] N.V. Krasnikov, A.A. Pivovarov, Phys. Lett. 112B (1982) 397 



17 



